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Abstract 

To understand the origin of hard X-ray emissions from black hole X-ray binaries during their low/hard 
states, we calculate the X-ray spectra of black-hole accretion flow for the following three configurations of 
hot and cool media: (a) an inner hot flow and a cool outer disk (inner hot flow model), (b) a cool disk 
sandwiched by disk coronac (disk corona model), and (c) the combination of those two (hybrid model). 
The basic features we require for successful models are (i) significant hard X-ray emission whose luminosity 
exceeds that of soft X-rays, (ii) high hard X-ray luminosities in the range of (0.4-30) xlO 37 ergs" 1 , and (in) 
the existence of two power- law components in the hard X-ray band with the photon indices of T s ~ 2 > Th, 
where T s and I\ are the photon indices of the softer (< 10 kcV) and the harder (> 10 keV) power-law 
components, respectively. Contribution by non-thermal electrons nor time-dependent evolution are not 
considered. We find that Models (a) and (b) can be ruled out, since the spectra are always dominated by 
the soft component, and since only one power-law component, at most, can be reproduced. Only Model 
(c) can account for sufficiently strong hard X-ray emissions, as well as the existence of the two power-law 
components, for a large ratio of the accretion rate in the corona to that in the thin disk. The outer disk 
corona (where the Compton y-parameter is smaller, y < 1) produces the softer power-law component with 
photon index of r s ~ 2, whereas the inner hot flow (where y £ 1) generates the harder component with 
Th < 2. This model can also account for the observed relationship between the photon index and the 
reflection fraction. 

Key words: accretion, accretion disks — black hole physics — hydrodynamics — radiation mechanisms: 
thermal 



1. Introduction 

It is well known that black hole X-ray binaries show 
two basic spectral states: the high/soft state appearing at 
high luminosities (typically 10% L^dd, where L^dd is the 
Eddington luminosity) and the low/hard spectral state 
at low luminosities, typically a few percent of £Edd, by 
factors of a few smaller than that of the high/soft state 
(McConnell et al. 2002). In the high/soft state X-ray 
spectrum is dominated by a thermal soft X-ray compo- 
nent of temperature, kT e ^ ~ 1 keV, associated with a weak 
steep power-law component with photon index T >2. In 
the low/hard state, on the other hand, spectra are mainly 
composed of Comptonizcd emission components in the 
hard X-ray band, soft excess components, and reflection 
components from the optically thick cold medium (for re- 
cent reviews, see McClintock & Remillard 2003; Done et 
al. 2007). While the high/soft state is reasonably well 
understood in terms of the standard disk model (see, e.g. 
Makishima et al. 1986), the nature of the accretion flow 
in the low/hard state is still a matter of debate. 

Historically, two classes of models have been consid- 
ered: the inner hot flow model and the disk corona model 
(Thornc & Price 1975; Shapiro et al. 1976; Liang & Price 
1977; for a recent review of theoretical models, see Kato 
et al. 2008). In the inner hot flow model, an optically 



thick disk is radially truncated at some radius, r tr , and 
the central region, r < r tr , is filled with a hot flow, which 
produces hard X-rays by means of the inverse Compton 
scattering of soft photons. The properties of the hot flow 
are modeled by an advection-dominated accretion flow 
(ADAF, or more generally, radiatively inefficient accretion 
flow, RIAF), in which most of the energy released via vis- 
cous dissipation remains in the accreting gas rather than 
being radiated away (Ichimaru 1977; Rees et al. 1982; 
Narayan & Yi 1994). The ADAF solution can exist only 
below a critical accretion rate, which depends on the vis- 
cosity parameter a (Ichimaru 1977; Narayan & Yi 1995; 
Abramowicz et al. 1995). Some authors calculated the 
radiative spectra from pure ADAFs or RIAFs by global 
Monte Carlo calculations (e.g. Kurpicwski & Jaroszynski 
1999; Kato et al. 2009). 

The disk truncation is suggested by a number of obser- 
vations. The soft excess component in the low/hard state 
is generally attributed to soft thermal emission from opti- 
cally thick, geometrically thin accretion disk (Balucinska- 
Church et al. 1995; Ebisawa et al. 1996). By spectral 
fitting of the soft thermal component to disk-blackbody 
(Mitsuda et al. 1984), the inner radius of the thin 
disk seems to be truncated at some outer radii, r tr ^ 10 
Schwarzschild radii (Makishima et al. 2008; Gierlihski et 
al. 2008). Moreover the low/hard state spectrum shows 
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a relatively small solid angle, 0, subtended by the reflec- 
tor derived from the reflection component (Gierlihski et 
al. 1997; Gilfanov et al. 1999). This implies that the in- 
terface between the hot flow and the cool disk has rather 
small area. One plausible explanation is that the inner 
part of the cool disk is truncated at some radius in the 
low/hard state. 

The disk truncation is also suggested by the variability 
study. The power density spectrum of the X-ray emis- 
sion shows a typical break frequency (Miyamoto et al. 
1992). As the spectrum evolves harder, the frequency 
becomes smaller (Ibragimov et al. 2005; Axelsson et al. 
2005). This would imply that the inner radius of the disk 
is far from the radius of the marginally stable orbit in the 
low/hard state, given that the typical frequency is related 
to the motion of the inner edge of the disk. 

In the disk corona model, by contrast, the optically 
thick disk remains untruncatcd, and hard X-rays are gen- 
erated in a hot corona. It is commonly assumed that a 
substantial fraction of the gravitational accreting energy is 
dissipated in the corona, although the physical processes 
by which the corona is heated are still unknown. Magnetic 
fields amplified inside the disk could dissipate above the 
disk surface, where the density is relatively low, making an 
active hot corona (Galeev et al. 1979; Haardt & Maraschi 
1991; Di Matteo 1998; Merloni 2003). The situation may 
be similar to the case of the solar corona. Liu et al. (2003) 
solved the radial structure of magnetically heated corona 
and calculated the emergent spectra by summing up the 
radiation from the plane-parallel corona at each radius. 
The corona could also be heated by the accretion of the 
coronal gas itself (Liu et al. 2002). Kawanaka et al. 
(2009) studied the structure of the hot flow by using non- 
radiative MHD simulation data, identifying it as a corona. 
They calculated self-consistent electron temperature and 
the emergent spectra by a global Monte Carlo calculation. 

More recently, however, it. is suggested by the Suzaku 
observations that the situation may be more complex. 
Takahashi et al. (2008) and Makishima et al. (2008) 
found that the hard X-ray continua of GRO J 1655-40 and 
Cyg X-l are better fitted by two thermal Comptonization 
components, rather than by single component. The exis- 
tence of the two Comptonized components in the hard X- 
ray band has already been previously reported (Gierlihski 
et al. 1997; Ibragimov et al. 2005). These would im- 
ply a possibility that a hot gas generating hard X-ray 
photons has a complicated structure with multiple values 
of electron temperature and/or optical depth. Another 
possibility of the origin of the two Comptonized compo- 
nents is a time-dependent effect (Poutanen & Fabian 1999; 
Takahashi et al. 2008). In this paper we consider the for- 
mer case and neglect the properties of the time variability. 

Though there are plenty of models for the hard X-ray 
spectra in the low/hard states, none of them can repro- 
duce the two power-law components so far. Especially 
the reproduction of the two power-law components is a 
theoretically challenging subject. In this paper we inves- 
tigate the structure and the emergent spectra of accretion 
flow for a variety of geometries of hot and cool accretion 



flows, instead of focusing on the detailed physical mecha- 
nisms making such accretion geometries, to examine which 
model can well explain the basic spectral properties of 
black holes during the low/hard state. Here we consider 
steady models and neglect non-thermal electrons which 
would explain the MeV tail in the low/hard state spectra. 
We will show that the combination of the inner hot flow 
model and disk corona model can naturally reproduce the 
two power-law components. In section 2 we present our 
models used for the spectrum calculations and details of 
the numerical method. Then we present the results of 
the emergent spectra in section 3. We discuss on the ob- 
tained results and applications for observation in section 
4. Finally, we give conclusions of this paper in section 5. 

2. Our Models 

2.1. Model Requirements 

In the present study, we aim at construction of accre- 
tion models which satisfy the basic spectral properties of 
the low/hard state. For this purpose we set the following 
three criteria for successful models for the low/hard state 
spectra (Done et al. 2007; Makishima et al. 2008): 

(i) The hard X-ray luminosity (at photon energy above 
3kcV), -Zj>3kcV7 should exceed the soft X-ray luminosity 
(below 3keV), L <3koV . 

(ii) The hard X-ray luminosity should be in the range 
of £> 3 kcv = 0.003 - 0.2L Ed d ^ 0.4 - 30 x 10 37 erg s" 1 for 
black hole mass of 10M Q . 

(iii) The spectrum should exhibit two power-law com- 
ponents which are dominant above and below ~ lOkeV, 
r s ~ 2 > Th, where T s (or Th) is the power-law index in 
the softer (harder) energy bands. 

Other prominent properties of the low/hard state, e.g. 
MeV tales in the spectra and time-lags in Fourier space, 
are not considered in this paper. 

2.2. Various Geometries of the Accretion Flows 

We consider three representative models for the config- 
uration of hot and cool flows (see figure 1). Generally 
speaking, hot and cool flows are expected to interact with 
each other by many different ways; e.g. energy exchange 
via radiation, conduction, and magnetic fields, mass ex- 
change via evaporation or condensation, and so on. 

Here, we focus on the radiative interaction of the flows 
and do not explicitly consider thermal conduction nor 
solve the dynamics of magnetic fields. Though there 
are some attempts for calculating the mass evapora- 
tion/condensation process between the hot flow and the 
cool thin disk (e.g., Rozahska & Czerny 2000; Spruit & 
Deufel 2002; Liu et al. 2007; Qiao & Liu 2009), we do not 
follow them since there exist uncertainties in the treat- 
ments, but instead incorporate the effects of matter ex- 
change by changing the accretion rate within the cool and 
hot flows. The exchanges of the energy and angular mo- 
mentum between the two phases (other than the energy 
exchange via radiation) are not considered in the present 
study. The effects of wind mass loss from the hot flow will 
be discussed in section 4. 
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The basic assumptions for each of the three models are 
as follows: 

Model (a): Inner Hot Flow Model 

In this model we postulate that a standard thin disk sur- 
rounds an inner hot ADAF. We assume that the material 
in the outer thin disk suddenly evaporates to the ADAF 
at the transition radius, r = r tr , so that the transition zone 
between the outer thin disk and the inner ADAF is suf- 
ficiently narrow. The accretion rate of the inner ADAF, 
Madaf, is assumed to be the same as that of the outer 
thin disk, Mdisk ; and these values are assumed to be radi- 
ally constant, Mtot = Madaf — Mdisk = constant, where 
Mtot is the total accretion rate of the hot flow and cool 
flow. Some fraction of the soft photons emitted from the 
thin disk is captured by the ADAF and cools the ADAF, 
especially near the transition radius, but only a small 
fraction of the hard photons generated by Compton up- 
scattering in the ADAF is caught by the outer thin disk 
because of the small cross section of the thin disk. The 
reprocessed radiation is thus expected to be weak. 
Model (b): Disk Corona Model 

In this model we consider two-temperature, accreting 
coronae above and below a thin cool disk. The corona 
is powered by the release of the gravitational energy by 
accretion itself. The thin disk extends down to the radius 
of the innermost stable circular orbit (ISCO) of the non- 
spinning black hole space-time, r = 3rs- We assume that 
the accretion rate of the corona, M coroIla , and that of the 
underneath thin disk, Mdisk, are both constant in space. 
Since the hot corona covers the whole thin disk, we expect 
that most of soft photons emitted from the thin disk expe- 
rience scattering in the coronae and that a large fraction 
of the hard photons generated in the coronae are repro- 
cessed by the thin disk. This makes a big contrast with 
the case of the inner hot flow model (a). 
Model (c): Hybrid Model 

We finally consider a hybrid model of Models (a) and 
(b). The outer part has a disk-corona structure, which 
is, however, truncated at some radius, r = r tr , and the 
flow becomes pure ADAF inside this radius. The mass 
accretion rate of the inner ADAF is a sum of that of the 
outer corona and the thin disk, i.e., Madaf = -Mtot = 
M corona + Mdisk- We assume that these values are con- 
stant in space unless specified otherwise. Note that similar 
configurations were considered in some previous studies, 
though the structure of the hot flow was assumed to be 
uniform (e.g. Poutanen ct al. 1997). 

2.3. Basic assumptions and simplifications 

In order to calculate the emergent spectra, we need to 
solve the radial structure of the ADAF, the corona, and 
the thin disk. Note that what matters in the thin disk 
part is only the radial profile of the soft photon emissiv- 
ity (or more precisely, the energy dissipation rate as a 
function of the radius) and it can be easily derived from 
the standard disk theory once an accretion rate is spec- 
ified. As for the structure of the ADAF and corona we 
adopt the technique used for the self-similar solution of 
the advection-dominated flow (Narayan & Yi 1995). Note, 
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Fig. 1. Schematic view of the three models of the accretion 
geometries. 

however, that the ADAF and the corona show somewhat 
different temperature profiles since the amount of soft ra- 
diation available for Compton cooling differs. 

In a realistic situation, dynamics of the magnetic fields 
and the way of interactions with the thin disk (via, e.g., 
thermal conduction and/or convective motions) could be 
markedly different in the corona and the ADAF. These 
issues are beyond the scope of the present study and are 
left as future issues. 

Throughout the present study we adopt the Newtonian 
dynamics and assume an axisymmctric and steady flow 
structure in the Newtonian gravitational potential, ip = 
—GM/r, where M is the black hole mass. We neglect 
the radiative force on the ADAF/corona, since the hot 
flow is optically thin and the typical luminosity is sub- 
Eddington. We further assume that the accreting gas is 
composed of two-temperature plasmas of fully ionized hy- 
drogen. Electron-positron pairs are not considered. 

Here the blackbody emission from the thin disk is as- 
sumed to be the source of soft photons and we neglect seed 
photons generated by synchrotron emission. Note that the 
synchrotron emission would be weak in a pure thermal 
plasma (Wardzinski & Zdziarski 2000). However the syn- 
chrotron emission would provide enough soft photons for 
the low/hard spectrum, if MeV(non-thermal) tail exists 
in the electron distribution, which is required in order to 
explain the observed spectrum in MeV band. (Wardzinski 
& Zdziarski 2001; Malzac & Belmont 2009; Poutanen & 
Vurm 2009). 

2.4. Basic Equations for Radial Structure of Hot 
Accretion Flow 

We employ the vertically integrated equations. 
Assuming the hydrostatic balance in the vertical direc- 
tion, we obtain the scale hight of the ADAF/corona, 



H 



vk 
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where c s = II/S is the sound speed, and E and II is verti- 
cally integrated gas density and pressure, respectively. We 
assume vertically uniform structure of the AD AF / corona 
within the half thickness, H . Hence the surface density 
is given by E = 2Hm p n, where n is the number density 
of protons, and the mass accretion rate of the ADAF(or 
corona) is M A DAF(M corona ) = -2nrT,v r . 

We obtain azimuthal velocity from the radial momen- 
tum equation, 

-— (2) 

Erflnr' [ ' 

where i>k is Kcplerian velocity and the last term is radial 
pressure gradient. In this paper we neglect radial inertial 
term for simplicity. It should be noted that this term con- 
tributes at most ~ 30% of the in equation (2) even for 
a = 1, which we adopt for reproducing the observational 
results. 

The angular momentum equation is 



2 2 



no radiation from inside this radius. In order to save 
the calculation time we restrict the computational region, 
r < Tout, where r out is the outer boundary radius and 
?"out = lOOrs unless stated otherwise. We employ loga- 
rithmically spaced 100 grid points to describe the prop- 
erties of the ADAF, the corona, and the thin disk. The 
basic equations can be integrated when cooling rates are 
calculated (see the next subsections). 

2.5. Coulomb Coupling and Bremsstrahlung Cooling 

The energy transport rate from protons to electrons by 
Coulomb collision is (Stepney & Guilbert 1983) 



d 1 d . 2 

E?v — (rvu,) = --r(r T rip ), 
dr r dr 



(3) 



where we assume that the r-ip component of the viscous 
stress tensor, T rip , is dominant and other components are 
negligible. We make the Ansatz that the stress tensor 
scales with the local gas pressure with a constant of pro- 



portionality, T rif 



-all. In order to solve the dynam- 



ical properties of the flow, we assume that the ADAF 
and the corona have self-similar structures, E ocr -1 / 2 and 

We assume that the electrons and protons are ther- 
modynamically coupled only through Coulomb coupling, 
neglecting poorly understood direct heating processes of 
electrons via magnetic fields. The energy balance between 
protons and electrons are, respectively, 

i du 7 ndE n d_ 

Coulomb " TJ-rip ^ 
' Qbrcms — Q Coulomb) \y) 



i - 



~l Vr ~ck 



7 . 



1 Vr E dr 



'Comp 



where the first two terms on left-hand-side of equation (4) 
represent advective cooling and the term on right-hand- 
side is viscous heating. Here an additional heating at 
the transition radius is not considered (Manmoto & Kato 
2000). We also assume that an advective term of electrons 
is negligible (see Nakamura et al. 1997). This means that 
the electron cooling timescale is much shorter than the 
accretion timescale, r/v r . We find that this assumption is 
valid unless the truncation radius is very large (> 100r s ). 
The vertically integrated pressure of the fully ionized 



gas is 



n = — 



(6) 



where T p is the temperature of protons and we neglect the 
pressure of electrons since their temperature is sufficiently 
low, T c <C T p . 

In our Newtonian model, the gas is assumed to be im- 
mediately swallowed by the black hole, once it reaches 
the inner boundary of the flow at 3rs, thereby emitting 
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where P and 9 C are the temperatures of the protons 
and electrons normalized by their rest mass energy, p = 
kTp/rripC 2 and O = kT c /m c c 2 . Here or is the cross sec- 
tion of the Thomson scattering, In A is the Coulomb log- 
arithm, In A ~ 20, and Ki is modified Bcsscl function of 
the second kind of order i. 

The bremsstrahlung cooling rate per unit area is 

Qbrcms = 2Hx afO- T m c c 3 n 2 F{O c ), (8) 

where «f is the fine-structure constant. Here the dimen- 
sionless factor is (Svensson 1984; Kato et al. 2008) 
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for 6 G < 1; 



for O > 1, 



C9) 



where 7e — 0.5772 is Eulcr's number. In our model we 
do not consider a supply of the soft photons from the 
bremsstrahlung emission in the Comptonization process, 
since the emission rate is sufficiently smaller than that of 
soft photons from the thin disk. 

2.6. Reprocessing of Comptonized Photons by the Thin 
Disk 

The dissipation rate of energy in the thin disk is 



/disk : 



3 gMMdisk 
47r r 3 



1 




(10) 



where r- ln = 3rs is the radius of the inner edge of the disk, 
where we adopt a torque-free condition. 

We take into account the radiative coupling between 
the corona and the thin disk by including reprocessing 
of radiation by the thin disk, as described by Haardt & 
Maraschi (1991). Some part of photons up-scattered in 
the ADAF/corona returns to the thin disk and heats up 
the surface of the disk. Hence the soft photon emission 
rate per unit area is the sum of the energy generated in 
the disk, Qdisk, and the reprocessed radiation, Q rep , 

(11) 
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where we set the albedo of the thin disk a = 0, i.e., all 
the photons impinging on the disk are assumed to be ab- 
sorbed. This soft photon emission cools down the elec- 
trons in the ADAF/corona by inverse Compton scatter- 
ing. 

2. 7. Monte Carlo Calculation 

We calculate the radiative spectra emerging from the 
accretion flow structures of the three geometries described 
above. We consider a computational region of a cylinder 
with a radius of r out and a height of 2r out . We neglect the 
height of the thin disk; we assumed that photons should 
be absorbed by the thin disk, if the path of the photon 
crosses the equatorial plane, z = 0, at r > rt r . We also 
neglect the absorption of the photon by free-free process. 
The bending of photon trajectory near the black hole nor 
absorption of photons by black hole event horizon are not 
included in this calculation. Inside the inner boundary of 
the hot flow, r < 3rg, we neglect any scattering. (In the 
computational procedure we set sufficiently large mean 
free path, A ~ 10 10 rs, and sufficiently low electron tem- 
perature, T ~ 10 eV.) 

We calculated the photon spectra scattered by hot 
electrons by using Monte Carlo calculation, following 
Pozdnyakov, Sobol, & Sunyaev (1977) and Liu, Mineshige, 
& Ohsuga (2003). We first choose the position of initial 
soft photon, r = ro, from the radial distribution of the soft 
photon cmissivity from the thin disk, Qsoft(^)- We set the 
weight wo — 1 for an initial photon with energy Eq, which 
has Planck distribution with temperature at ro. We cal- 
culate the probability Pq = cxp(— r) for photon passing 
through the hot gas, where r is optical depth evaluated 
from the photon position to the boundary of computa- 
tional region or the thin disk. Then wqPq is the transmit- 
ted portion of photon and the remaining w\ = wq(1 — Pq) 
is the portion of photon scattered at least once. Let 
w n = tu n _i(l — P n -i) be the portion of photon experi- 
encing the n-th scattering. We continue the calculation 
until the weight becomes sufficiently small, w n < 10~ 7 . 
Repeating the same procedures for sufficiently large num- 
ber of photons, typically I0 8 , we calculate emergent spec- 
tra by collecting photons escaping from the computational 
region and reprocessed radiation by collecting photons ab- 
sorbed by the thin disk. We also calculate the cooling rate 
of the electrons at each radial bin by evaluating the change 
of the photon energy at every scattering. In this paper we 
do not calculate the reflection features from the thin disk. 

By this Monte Carlo calculation, we obtain the quan- 
tities of radiative processes, Qcomp and Q r cp- In order 
to obtain self-consistent structures of accretion flows, we 
use these Qcomp and Q rcp values to re-calculate the flow 
structure and iterate (see below). 

2.8. Iterative Method 

We first solve the structure of the ADAF and the 
corona by using tentative values of the Compton cooling 
rate, Qcomp.cb and the emissivity of reprocessed photons, 
Qrcp.Oj where the subscripts denote the number of itera- 
tions. Then we solve the radiative transfer by the Monte 
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Carlo calculation, obtaining the exact values, Qcomp,i 
and Q rep .i, for this structure. We solve the structure 
of the flow again by using the newly obtained cooling 
rate, Qcomp.i, and calculate the radiative spectra by using 
Q rC p.i. Repeating the same procedure, we finally obtain 
the self-consistent structure and radiative spectra. 

3. Results of the X-ray Spectra 

There are four main model parameters: the viscos- 
ity parameter a, the total mass accretion rate mtot, 
TOcorona/^diskj and the truncation radius r tr , where rhtot = 

M tot C 2 1 XEdd, "Icorona = ^coron a C 2 /iEdd, and TOdisk = 

-MdiskC 2 /-/jEdd- For clarity we show the typical values of 
parameters used in this paper: a =0.3 and 1, rhtot =0.1, 
0.3, 0.8, and 1, m corona /7« d isk = 1 - 1000, and r tr /r s = 
3 — 300. Note that relatively large values of a is neces- 
sary to obtain the advection-dominated solutions at large 
accretion rates (i.e., hard X-ray luminosities). The black 
hole mass, m = M/M Q , should be adjusted to each ob- 
ject. Hereafter we study the cases of the stellar mass black 
holes, m = 10. 

3.1. (a) Inner Hot Flow Model 

In this model the mass accretion rate of the ADAF, 
mADAF = A^ADAFC 2 /iEddj is assumed to be the same as 
that of the thin disk. Note that jtiadaf cannot be larger 
than the critical value, which depends on the viscosity 
parameter a. 

We calculate the radial structure of the ADAF and show 
the results for a = 1, jtiadaf = 1, and r tr = 30rg in fig- 
ure 2a. Here the vertical optical depth is th = na^H, 
and the Compton y-parameter is y = (48 G + 168 2 )t# . 
The advection fraction, / a d v = Qadv/Qvisc, is approxi- 
mately unity, validating the assumption that the hot flow 
is advection-dominated. This means that the Coulomb 
coupling timescale is much longer than the accretion 
timescale. Hence the scale height and the optical depth 
follows the advection-dominated solutions, H ~ r and 
th oc a~ 1 mADAF?" _1 ' 2 as long as / a d v ~ 1 (Narayan & 
Yi 1995). The ratio of the bremsstrahlung emission rate 
to the Compton cooling rate, / bc = Qbrems/Qcomp, is so 
small that we can safely neglect the contribution of the 
bremsstrahlung emission to the spectrum. 

We display the emergent spectra from this structure in 
figure 2b. Here we do not consider the reflection feature in 
the spectrum, which is not expected to be so strong. Since 
the Compton cooling is most efficient at r = r tr , where the 
emission of the soft photons from the thin disk is largest, 
the electron temperature decreases with radius. Moreover 
the optical depth increases with decreasing radius, so does 
the y-parameter. Since the soft photons come only from 
the thin disk part located at r > r tr , a small fraction of the 
soft photons is intercepted by the inner ADAF, resulting 
in a very weak Compton component in the radiative spec- 
tra (about 30% of the thermal component; see Table 1), 
which does not agree with the observations. Note that we 
only consider the soft photons from the thin disk inside 
the boundary r < r out = 100rs. Thus we slightly undcrcs- 
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Fig. 2. The structure of the ADAF (2a: left panel) and the emergent spectra (2b; right panel) for a = 1, ritADAF = 1) an d rtr = 30rg . 
Left panel (2a): In the ADAF region, 3rg < r < rtr, we show the dimensionlcss values, electron temperature, O e (solid line), vertical 
optical depth, rjj (long-dashed line), Compton {/-parameter (short-dashed line), advection fraction, / a( j v (dash-dotted line), and the 
ratio of the cooling rate of brcmsstrahlung emission to the Compton cooling rate, /t> c (dotted line). Right panel (2b): The spectrum 
of the photons escaping from the computational region and soft photons from the thin disk is shown by solid and dashed line, 
respectively. 



timate the soft X-ray luminosity since we neglect the soft 
photons from the thin disk at r > r out (see below). 

We calculate variety of the emerging spectra for other 
parameters (Table 1). We find that the hard X-ray com- 
ponent is always much weaker than the soft X-ray compo- 
nent. As rtr decreases, the ratio of the hard X-ray lumi- 
nosity to the soft X-ray luminosity, £>3kev/i<3keV, de- 
creases, though the most of the energy is dissipated in the 
inner part of the ADAF. This is because an increase in 
the number of soft photons results in decreasing electron 
temperature due to more efficient Compton cooling of the 
ADAF region. 

We also calculate the radiative spectra for much larger 
values of the truncation radius than that estimated by the 
observational results, r tr ~ 10rg (see e.g. Makishima et al. 
2008). In such cases we should set large outer boundary 
r out = lOOOrs in order not to underestimate the soft pho- 
ton luminosity from the outer thin disk. We see that the 
radio of the hard X-ray luminosity to the soft X-ray lu- 
minosity is still small (< 0.3) for large values of r tr (Table 
1). To make the situation worse, the luminosity of hard 
X-ray, as well as soft X-ray, decreases as the truncation 
radius increases. This is because, the larger r tr is, the 
smaller becomes the number of the soft photons entering 
the inner region of the ADAF and the less efficient be- 
comes Compton cooling. Then the electron temperature 
becomes higher and the energy transfer rate from pro- 
tons to electrons by Coulomb collisions (which is approxi- 
mately proportional to To 3 ^ 2 ) becomes smaller (Coulomb 
coupling timcscalc becomes longer). Consequently the 
hard X-ray luminosity decreases (advection fraction in- 
creases) as rt r increases. The hard X-ray luminosity in 
the case of r tr = 300rs is one order of magnitude smaller 



than the criterion for the low/hard state. Note that for 
large r tr (^ 300) the bremsstrahlung emission from the in- 
ner region of the ADAF dominates over the Comptonized 
radiation. To summarize, we cannot reproduce the X-ray 
spectra in the low/hard state by the inner hot flow model 
for (a,m A DAF) = (1,1). 

For a given a, the hard X-ray component becomes 
harder and stronger as the accretion rate, rh-ADAF in- 
creases. However, the mass accretion rate cannot be larger 
than the critical value, rh cr i t ~ l(or 0.1) for a = 1(0.3) in 
the inner ADAF. We can see that even if the accretion 
rate is its maximum value for a given a, the hard X-ray 
component is still weaker than the soft X-ray component 
(see the case of a = 0.3 in Table 1). 

In summary the inner hot flow model always produces 
too weak hard X-rays to explain the observed spectra in 
the low/ hard state. 

3.2. (b) Disk Corona Model 

In this model we assign three parameters: the viscosity 
parameter a, total accretion rate rh, to t = "^corona + ^disk, 
and the ratio of the accretion rate of the corona to that of 
the thin disk rh, cor0 na/TOdisk- We show the structure of the 
coronal flow for (a, m to t, ra corona /mdi s k)=(l, 0.8, 79) in 
figure 3a. The maximum accretion rate for a = 1 is found 
to be m t ot ~ 0.8 in the disk corona model. The electron 
temperature is almost constant in space, in contrast to 
the case of the inner hot flow model. This is because 
significant Compton cooling in the innermost region by 
the underlying thin disk where cmissivity rapidly increases 
inward. Nevertheless the inner region has the largest y- 
parameter, because of a large optical depth. Note that the 
electron temperature increases at the inner edge, since the 



No. ] Radiative Spectra in Black Hole X-ray Binaries 7 



Table 1. Parameters and the resultant values for the inner hot flow model 



a 


"1ADAF 


n T /r$ 


£<3kcv [erg s" 1 ]* 


£>3kcv [erg s -1 ]* 


ratio* 


r* 


1 


1 


10 


4.2 x 10 aY 


7.2 x 10 3b 


0.17 


2.05 


1 


1 


30 


1.4 x 10 37 


4.6 x 10 36 


0.33 


1.88 


1 


1 


30* 


2.0 x 10 37 


4.8 x 10 36 


0.23 


1.88 


1 


1 


100* 


6.2 x 10 36 


2.0 x 10 36 


0.32 


1.88 


1 


1 


300* 


1.7 x 10 36 


6.1 x 10 35 


0.35 


1.70 


1 


0.3 


30 


3.9 x 10 36 


2.3 x 10 35 


0.06 


2.11 


0.3 


0.1 


30 


1.3 x 10 36 


2.3 x 10 35 


0.17 


1.98 



* The integrated luminosity below and over 3 keV, and the ratio £>3k c v/^<3kcV- 
t Photon index evaluated by two points of 3 and 100 keV in the spectra. 

* r Du t = lOOOrs- 



emission rate of the soft photons, g so ft = r 2 Q so { t / Mjot^Ki 
vanishes at the torque free radius. 

The solid line in figure 3b shows the result of the spec- 
trum from the coronal structure described above. The 
emergent spectrum has a single power-law component in 
the hard X-ray band, though the y parameter varies in the 
radial direction. We find that the hard X-ray luminosity 
is only ~ 40% of the soft X-ray luminosity, and the photon 
index is T ~ 2.2. Hence the spectrum does not satisfy the 
criteria for the low/hard state, even when a large ratio of 

TO corona /TOdisk IS adopted. 

For small mass accretion rate, rh tot = 0.3, the resulting 
spectrum has two humps in the hard X-ray band, rather 
than a single power-law component (see figure 3b). This 
is because the electron temperature is so high, Q e ~ 1, 
that the photons gain significant energies only by a sin- 
gle scattering. We can see that the hump around 10 keV 
is made by one scattering of the soft photons, and an- 
other high energy hump is made by the photons scattered 
at least twice. We also show the emergent spectra for 
(a, rhtot) = (0.3,0.1) in figure 3b. We find that its spectral 
features are very similar to those for (a,rh to t) = (1,0.3). 
Note that we cannot find the advection-dominated solu- 
tion of the corona with larger accretion rate, rhtot > 0.1, 
for a = 0.3. 

We show in figure 4 (squares) the resultant luminosi- 
ties and photon indices for various ratios of the accretion 
rate rh CO rona/rhdisk but keeping rhtot = 1- We see that the 
spectrum has harder (smaller T) and stronger power-law 
component for larger ratios of rh CO rona/rhdisk. However 
the hard X-ray luminosity is less than 50% of the soft X- 
ray luminosity and the photon index is larger than 2, even 
when rhcorona/rhdisk = 1000. 

The situation is not improved but rather getting worse 
for the cases with a smaller rh to t = 0.3 or with a smaller 
a = 0.3 (see triangles and diamonds in figure 4). In the 
both cases, the hard X-ray luminosity is much smaller 
than the observed values, L>3kcV ^ 2 x 10 35 erg s _1 . 

For the coronal geometry we expect that the fraction of 
the escaping photons from the corona should be compa- 
rable to that of the reprocessed photons. Hence, the ther- 
mal component should exceed the power- law component. 
Similar results were obtained previously by an one-zone 
disk corona model (e.g. Haardt & Maraschi 1993; Stern 
et al. 1995; Poutanen et al. 1997; Dove et al. 1997; Cao 



2009). We confirm these results by the radially structured 
corona models. We conclude that the hard X-ray compo- 
nent can never exceed the soft X-ray component even for 
large rh to t and large rh cor0 na/rhdisk and the photon index 
of the power-law component is always lager than two. The 
situations will be improved, if the thin disk is truncated 
(see below). 

3.3. (c) Hybrid Model 

We study the hybrid case, in which a thin disk covered 
with disk coronae is truncated at r = rtr and transformed 
to the ADAF . The structure of the ADAF/corona is 
shown in figure 5a for (a,mtot,rh corona /rhdi s k,r'tr/r's)=(l, 
1, 99, 10). We can see that the radial gradient dis- 
continuously decreases inside r tI where the underneath 
thin disk does not exist. Because of the differences in 
the radial distribution of the electron temperature (and 
thus the y-parameter), the outer corona and the inner 
ADAF produce hard X-ray with different power-law in- 
dices. The spectrum has a relatively flat power-law com- 
ponent (r s ~ 2) in the intermediate energy range, 1 — 10 
keV, generated by the outer corona with small y < 1, and 
a hard power-law component (T^ < 2) in the high energy 
range, > 10 keV, with cut-off at a few hundreds keV, by 
the inner ADAF with y > 1 (figure 5b). Unlike the in- 
ner hot flow model, the hard X-ray component can be 
dominant, as long as the accretion rate of the corona is 
sufficiently larger than that of the thin disk (typically 
rh C orona/rhdisk ^ 100), since the thin disk emits less soft 
photons than that in the inner hot flow model with the 
same total accretion rate. 

We show the features of the emergent spectra for vari- 
ous parameters in figure 6. We cannot find the advection- 
dominated solutions for small values of rt r with rhtot = 1 
because of efficient Compton cooling. For r tr ^ 6rs we can 
see the trend that the larger rt r is, the smaller becomes 
the photon indices of the power-law components and the 
larger becomes the difference between the two photon in- 
dices (r s — Th). This means that the break between the 
two power-law components is more prominent, when the 
thin disk is truncated at a larger radius. 

We see that the soft X-ray luminosity steeply decreases 
as rt r increases, while the hard X-ray luminosity only 
gradually decreases. Thus the hard X-ray component 
can dominate over the soft one for large r tr (> 10rg for 
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Fig. 3. The structure of the corona (3a; left panel) and the emergent spectra (3b; right panel). Left panel (3a): We show 
C (solid line), th (long-dashed line), and y(short-dashed line), as well as the soft photon emissivity, q so ft (dash-dotted line) for 
(o,mtot,m corona /m disk )=(l, 0.8, 79). Right panel (3b): Emergent spectrum for (a,m t ot,m corona /m disk )=(l, 0.8, 79) (solid line), 
(1, 0.3, 299) (long-dashed line), and (0.3, 0.1, 99) (short-dashed line). 




Fig. 4. The spectral features of the disk corona model as functions of the ratio of the accretion rate of the corona to that of the 
thin disk, m CO rona/"idisk> f° r given a and rhtot- Squares, triangles, and diamonds represent the results for (a,rhtot) = (1,0.8), 
(a,mtot) = (1,0.3), and (a,mtot) = (0.3,0.1) respectively. Left panel (4a): The luminosity of the emergent spectra below and over 
3keV, £<3kcV (open symbols) and L>3keV (filled symbols). The spectrum satisfies the criterion for -L>3k e v when the filled symbol 
is above the shaded region. Right panel (4b): The photon index evaluated by the two points, 3keV and lOOkeV. The photon index 
should be below the shaded region in order to satisfy the criterion. 
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Fig. 5. The structure of the accretion flows (5a; left panel) and emergent spectrum (5b; right panel) for (a,mtot ,^corona/"idisk) = (l; 
1, 99) in the hybrid model. The thin disk is truncated at rt r = 10rg. Left panel (5a): Same as figure 3a but for the hybrid model. 
Right panel (5a): Emergent spectrum (solid line) and soft photon spectrum from the thin disk (dashed line). 




Fig. 6. The spectral features of the hybrid model for a = 1 and rritot = 1 as functions of the truncation radius, rt r . The ratio of 
the accretion rates are m CO rona/™.disk = 1 (diamonds), m CO rona/«idisk = 9 (triangles), and m coron a/mdisk = 99 (squares). Left panel 
(6a): The luminosity L>3keV an d ^<3keV> i s plotted by the filled and open symbols, respectively. The hard X-ray luminosity should 
be above the shaded region for the low/hard state spectrum. Right panel (6b): Photon index of the softer one of the power-law 
components, r s (evaluated by the two points of 3keV and lOkeV), and that of the harder one, Th (evaluated by the points of 30kcV 
and lOOkeV). The spectrum satisfies the criterion if r\ is below the shaded region. 
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TOcorona/^disk = 99). This is because the reprocessed 
radiation becomes weaker when the disk is truncated 
at a larger radus. For a very large r tr (3> 10rs for 
Wcorona/Wdisk = 99), on the other hand, the absolute value 
of the hard X-ray luminosity becomes less than that of the 
criterion. The ratio of the luminosity £>3kev7-k<3keV de- 
creases as the ratio of the accretion rate m corona /mdisk 
decreases, and for small m corona /TOdisk < 10, it is hard to 
obtain the results satisfying the criteria for the low/hard 
state. For a small total mass accretion rate rhtotj the hard 
X-ray luminosity is significantly reduced as we saw in the 
disk corona model (section 3.2). Finally, we find that the 
resultant spectra of the hybrid model agree well with the 
observational results only for m to t = 1, "icorona/wdisk^ 10, 
and r tI ~ 10 — 20rs, though the typical electron temper- 
ature is 200-300 keV, which is slightly larger than that 
estimated by the X-ray observation, ~ 100 keV. 

4. Discussions 

4-1- General Discussion 

4- 1.1. Brief Summary 

We calculated the expected spectra from the black hole 
accretion flows with three different flow geometries. We 
find that only the hybrid model can account for the ob- 
served spectra in the low/hard state for the following rea- 
sons, (i) It can reproduce the two power-law components 
from the two separate hot regions having different ther- 
mal properties: the outer corona and the inner ADAF. 
The former produces the softer power-law component with 
r s ~ 2, while the latter generates the harder one with 
Ft < 2. (ii) The spectrum is dominated by hard X-ray ra- 
diation because of less soft photons originating from the 
outer thin disk impinging the two hot zone (i.e., we find 
less dissipated energy in the thin disk compared to the 
inner hot flow case and less reprocessed radiation com- 
pared to the disk corona case). When the number of the 
soft photons is small, the soft thermal emission becomes 
weak and the Comptonization in the hot zone becomes 
very efficient since y is moderately large (y > 1). 
4-1.2. On the a Value 

When the optical depth is much smaller than unity, 
for a moderately large y ~ 1, the electron temperature 
should be very high, kT e > m c c 2 . However, the resulting 
spectrum does not show a simple power-law component 
with such a high electron temperature (see figure 3b). 

Therefore the spectrum can have a strong power-law 
component only for a moderately large optical depth, 
Tff > 1. Since the optical depth of the hot advection- 
dominatcd flow is given by r# -~ m/a, the smaller a is, the 
larger becomes the optical depth for a given accretion rate. 
However the solution of advection-dominated flow exists 
only when the mass accretion rate is smaller than critical 
value, m cr ;t oc a 2 (Narayan et al. 1995; Abramowicz et 
al. 1995). Hence the maximum optical depth is larger 
for larger a, th ct, and large a(^ 1) is necessary for a 
moderately large optical depth of the hot flow, th ~ 1, to 
be realized (Zdziarski 1998). 

Global three-dimensional MHD simulations, on the 



other hand, showed that the effective value of the vis- 
cosity parameter is very small, a = 0.02 — 0.2, in the non- 
radiative ADAF or RIAF (Hawley et al. 2001; Machida 
& Matsumoto 2003). However the ADAF/corona which 
we consider is radiative and interacts with the thin disk. 
Local radiative MHD simulation shows magnetically dom- 
inated corona forming above the thin disk (Miller & Stone 
2000; Hirose et al. 2006). Moreover, when the hot flow 
collapses in the vertical direction by radiative cooling, 
we could expect a magnetically supported accretion flow 
(Machida et al. 2006; Oda et al. 2009). In the flow sup- 
ported by the magnetic fields, the effective a could be 
large, a ~ 1, though our model does not take into account 
the magnetic processes explicitly, and the a prescription 
might not be relevant in such a low /3(< 1) region. 
4-1.3. Heating of ADAF/corona by Magnetic Fields 

Since the hard X-ray photons are emitted from the 
ADAF/corona regions where the advection fraction is al- 
most unity, the radiative efficiency of the flow we consider 
is very small, L/M tot c 2 ~ 0.01. This conclusion would be 
changed, if there are direct heating of electrons by mag- 
netic processes, in addition to the electron heating by the 
Coulomb collisions with protons. Moreover since the heat- 
ing by magnetic fields is non-local process, the energy dis- 
sipated in some region can be transported to other distant 
region via magnetic fields. 

In the inner hot flow case, if a significant fraction of en- 
ergy released in the outer thin disk is dissipated in the in- 
ner ADAF region, the radiative efficiency of the thin disk 
can be sufficiently small (< 1), and that of the ADAF 
becomes large (< 1). As a result, the hard X-ray luminos- 
ity from the ADAF can exceed the soft X-ray luminosity 
from the thin disk. Dove et al. (1997) and Poutanen et 
al. (1997) considered similar situation and claimed that 
inner hot flow model can basically explain the low/hard 
state spectra. 

In the disk corona case, on the other hand, the con- 
clusion that soft X-ray luminosity is always dominant 
is hardly affected by the heating process of the corona. 
Even if all of the accretion energy of the disk and corona 
is released by the hard X-ray radiation in the corona, 
half of the hard X-rays is reprocessed by the thin disk. 
Consequently the hard X-ray luminosity cannot be larger 
than the soft X-ray luminosity (Dove et al. 1997; Cao 
2009). 

4-1-4- Mass Loss from the ADAF/corona 

So far we assumed no mass loss, thus setting m to t=const 
in the radial direction. We next calculate the effect of the 
mass loss from the ADAF/corona by using radially depen- 
dent mass accretion rate in the case of the hybrid model. 
For simplicity, we assumed that the total accretion rate 
decreases with decreasing radius, m to t = "T-tot.out^Aout)* 
with s being a constant, while m corona /rhdi s k is kept con- 
stant. Here m to t,out is the total mass accretion rate at the 
outer boundary, r = r out (= 100rs). In figure 7 we show 
the results of the emergent spectra for various values of s. 
We find that as s increases the total luminosity decreases 
and the power-law component becomes softer and weaker 
compared with the thermal component. Thus the effect 
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Fig. 7. The emergent spectrum of the hybrid model for 
(a, rh corona /rh disk , r tr /r s )=(l, 99, 10). The total mass 
accretion rate is assumed to be rhtot = "itot,out(f/f'out) s 
where mtot.out = 1 and s=0(solid line), 0.1 (long-dashed line), 
0.5 (short-dashed line), and 0.8(dash-dotted line). 

of the mass loss only weakens the hard X-ray component. 
If only the accretion rate of the corona decreases with 
decreasing radius while that of the thin disk is constant 
(TO cororia /rhdisk decreases), we can expect that the hard 
X-ray component becomes even weaker than that in the 
above case. Therefore the effect of the mass loss results 
in an unfavorable consequence for the low/hard state. 

4-2. Observational Implications 

4-2.1. Photon Index and Truncation Radius 

Zdziarski ct al. (1999) found the correlation between 
the indices of the power-law component and the reflection 
fractions, R, in the accreting systems. They explained 
this correlation by changing the inner radius of the outer 
thin disk surrounding an inner hot flow. Let us confirm 
their result by our hybrid model. 

Since we do not include the reflection features in the 
spectra and it is difficult to estimate the reflection frac- 
tions in the complicated accretion structures as are en- 
visaged by our models, in order to calculate the reflec- 
tion fraction for a given r tr , we thus consider a simple 
model. Let us consider a point source of hard X-rays lo- 
cated above the central black hole, (r, z) = (0,h), where h 
is the heigh of the hard X-ray source. We calculate the 
solid angle, f2, subtended by the thin disk which is trun- 
cated at rt r and which extends to infinity at the outer 
part. Then we obtain the reflection fraction, R = f2/27r. 
In figure 8 we show the relationship between R and Fh 
calculate based on our model, overploting the observed 
values of Cyg X-l and GX 339-4 taken from Zdziarski et 
al. (1999). We find that our model reasonably agrees with 
the observed correlations, especially for h = 5rg, though 
the data are quite dispersed. Note that even if the disk 
is truncated at the radius far from the black hole, signifi- 
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Fig. 8. Correlation between the photon index I\ 
and reflection fraction R in the hybrid model for 
(a, m t ot, m C orona/»Ttdisk)=(l! b ln order to calcu- 

late R, we assume the height of the hard X-ray source, 
h = 10rg (squares) , h = 5rg (triangles) , and h = 3rg (diamonds) . 
The data points are taken from Zdziarski et al. (1999). 

cant fraction of the hard X-rays originates from the inner 
ADAF region. Hence we expect that the typical position 
of the hard X-ray does not change so significantly, even 
when r tr increases. 
4-2.2. Two Power-law Components 

We can reproduce the two power-law components, 
which arc suggested by the observations. The softer 
power-law originates from the outer disk corona, where 
y is relatively small, while the harder one originates from 
the inner ADAF. The photon indices of the two power- 
law components obtained in our model are (r s ,Fh) ~ 
(2.0,1.8) for (a,mtot,"icorona/TO disk ,rt r /rs)=(l, 1, 99, 10), 
which roughly agree with the observed values of Cyg X-l, 
(r s ,r h )~ (1.9,1.7), and GRO J 1655-40, (r s ,r h )~ (1.8,1.6) 
(Takahashi et al. 2008; Makishima et al. 2008). 

Note that the two power-law components are produced 
given that the thin disk is truncated at some radii and 
that the ADAF/corona are interacting with the soft pho- 
tons from the thin disk. In other words the appearance 
of the two power-law components is a good indication of 
the disk truncation. We see that the spectral features de- 
pend on rt r (see figure 6), so the spectral evolution of the 
power-law component provides information regarding the 
evolution of the truncation radius; the larger r tr is, the 
larger become the difference between the photon indices, 

r s -r h . 

Note that there is another explanation of the origin of 
the two power-law components. It is well known that the 
spectrum varies on short time scales. It is possible that 
two power-law components could be artificially generated 
because of the time average in the long observational time. 
If one tries to explains also time lags, then one would have 
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to conclude that either the local spectrum is pivoting due 
to evolution of the geometry of the energy dissipation re- 
gion (Poutancn & Fabian 1999), or the spectrum changing 
because at different distances from the black hole the spec- 
trum of the hot flow is different and it is harder closer in 
(Kotov et al. 2001). 
4-2.3. Soft Thermal Component 

We evaluate the soft thermal component from the thin 
disk, including the reprocessed radiation. It should be 
noted that the spectrum of the soft photons in our model 
is significantly different from the disk-blackbody spectrum 
with an accretion rate of Mdisk- Rather it exhibits a 
high effective temperature (or large spectral hardening 
factor, k), because the contribution of the reprocessed 
radiation of the hard X-ray from the ADAF/corona is 
large (typically ^ 50% of the soft photons). Hence, the 
simple estimation of the inner disk radius by the disk- 
blackbody model would underestimate the actual disk ra- 
dius (Gierlihski et al. 2008). 
4-2.4- Spectral Transitions 

We find that the spectra show very different spectral 
features for various ratios of the accretion rate of the 
corona to that of the thin disk, even when the total ac- 
cretion rate is kept constant. The larger TOcorona/^disk is, 
the larger becomes the fraction of hard X-rays (see figure 
6a). This implies that a change in the ratio of the ac- 
cretion rate may trigger a transition of the spectral state, 
even when the total mass accretion rate is nearly kept 
constant. Note that our model produce a thermal cut- 
off in the power-law component around ~ 100 keV. Thus 
the non-thermal electrons would be necessary for the spec- 
trum in the high/soft state, and also in the low/hard state 
if the observed radiation in the MeV band comes from the 
accretion flow (Li et al. 1996; Poutanen & Coppi 1998; 
Gierlinski et al. 1999; McConnell et al. 2002; Yuan et al. 
2003; Inoue et al. 2008). 

4-3. Further Developments of Our Model 

The electron temperature expected by our model (typ- 
ically 200-300keV) is somewhat higher than that esti- 
mated by the cut-off energy of the hard X-ray spectrum 
(~100keV). This is because, even though the viscosity pa- 
rameter is large, a ~ 1, and the mass accretion rate is 
its maximum value m to t, the optical depth of the ADAF 
is still somewhat smaller than unity, resulting in high 
electron temperature in order to keep large Compton y- 
paramcter, y ~ 1 (Zdziarski 1998). The optical depth of 
the electron-positron pairs could contribute to lower the 
electron temperature, since the compactness parameter 
around the black hole in the low/hard state is expected 
to be large, 

, S -^L~10(— ^- )(^-Y\ (12) 

where we set the typical size of the emitting region, 
R = 10rg. The newly found solutions, luminous hot ac- 
cretion flows, might also be appropriate, since they have 
large optical depth and low electron temperatures, com- 
pared with ADAFs for a given y (Yuan 2003; Yuan k, 



Zdziarski 2004), though they might be thermally unsta- 
ble. Moreover thermal synchrotron emission would play a 
role as a cooling mechanism of the electrons, in addition 
to a role as a source of the soft photons. 

We find that the hybrid model with the parameters used 
in this paper can marginally agree with the relatively faint 
low/hard state (L/LEdd ~ 0.01). However the hard X-ray 
luminosity obtained by our model is still fainter than the 
relatively luminous low/hard state, L/ivEdd^0.1. In order 
to explain this state, we could consider other accretion so- 
lutions of the ADAF/corona: e.g., luminous hot accretion 
flow (Yuan & Zdziarski 2004) or magnetically dominated 
accretion flow (Oda et al. 2009). 

In our model we simply give the transition radius and 
mass accretion rate as free parameters. However these 
values seem to depend on the Compton cooling rate in 
the disk evaporation model (Liu et al. 2002). Hence we 
have to solve the whole of the structure, taking in to ac- 
count the gas evaporation/condensation process and ra- 
diative transfer self-consistently. In addition we used the 
self-similar solution for ADAF/corona, which would be 
invalid near the black holes where boundary conditions 
become essential (Nakamura et al. 1996; Manmoto et al. 
1997; Narayan et al. 1997). For further studies, we should 
include general relativistic effect into both accretion struc- 
ture and radiative transfer (doppler beaming by the gas 
motion and bending of the photon trajectory), in order 
to include the radiative contribution from the plunging 
region where the bulk Comptonization would also be im- 
portant. 

5. Conclusions 

We calculate the emergent spectra from three accre- 
tion geometries: (a) inner hot flow model, (b) disk corona 
model, and (c) hybrid model. We find following results. 

1. Only the hybrid model satisfies all the basic features 
of the low/hard state for m corona ~ 1, m disk < m C orona, 
and rt r ~ 10rs. Other models predict weak hard X-ray 
emission and at most one power-law component. 

2. Two power-law components simply imply that the 
thin disk is truncated and there are two regions with differ- 
ent thermal properties: the outer corona produces softer 
component dominating at < 10 keV and the inner ADAF 
generates a harder component. 

3. The hybrid model can also explain the observed T — R 
relation. 
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